Evidence for ideal insulating/conducting state in a ID integrable system 
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Using numerical diagonalization techniques we analyze the finite temperature/frequency conduc- 
• tance of a one dimensional model of interacting spinless fermions. Depending on the interaction, 

' the observed finite temperature charge stiffness and low frequency conductance indicate a funda- 

, mental difference between integrable and non-integrable cases. The integrable systems behave as 

ideal conductors in the metallic regime and as ideal insulators in the insulating one. The non- 
integrable systems are, as expected, generic conductors in the metallic regime and activated ones in 
the insulating regime. 
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PACS numbers: 05.45. +b, 71.27.+a, 72.10.-d 



In classical integrable systems the existence of a macroscopic number of conservation laws has a profound conse- 
quence, dissipationless transport 0]. The analogous effect in nontrivial quantum integrable systems has been studied 
only recently in a prototype model of dissipation M . In this model, describing a particle interacting with a fermionic 
bath in one dimension, we found that the tagged particle shows ideal mobility even at finite temperatures, T > 0, 
when the system is integrable. It is desirable to extend these ideas to homogeneous many-body models, which are at 
present of particular interest in connection with strongly correlated electrons. In this field, most analytical findings 
are on integrable ID models as the Hubbard || or the spinless fermions with nearest neighbor interaction model [|J. 
Progress in the study of dynamical response at T > is hindered by the lack of reliable methods. The only attempts, 
' which however might obscure the role of integrability, start from a Luttinger liquid effective Hamiltonian description 

^ ' In this work, based on a recent reformulation of the finite temperature charge stiffness El and numerical methods 
for calculating T > dynamical conductivities cr(u>) on finite size systems, we present evidence of the importance 
of integrability for transport properties. 

From linear response theory, the real part of the conductance cr(u>) at frequency to is given by: 



<j(u}) = 2ttDS(u>) + a reg (ui), (1) 

1 — e _/3w 7r . 

OYe 9 M = y V] p n \ (n \ j \ m) \ 2 6(lu - e m + e n ), (2) 
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where | n),e n denote the eigenstates and eigenvalues of the Hamiltonian, p n the corresponding Boltzmann weights, j 
the current operator and f3 = 1/fceT. We will consider ID tight-binding models of L sites (ks = h = e = 1). 
cr{ui) satisfies the optical sum rule 0: 

a(uj)dw = -(-f), (3) 

oo L 

where (T) denotes the thermal expectation value of the kinetic energy. The sum rule (|) together with Eqs.(|J|) can 
be used for the evaluation of the stiffness D(T) ^|J2|. It will be however more convenient to discuss the behavior of 
D at finite temperatures, with a generalization of the original Kohn's approach Q for zero temperature, by relating 
D(T) to the thermal average of curvatures of energy levels subject to a fictitious flux <j) j| : 

n 1 n 1 d 2 e n {4>) , , . 

D = T l^PnD n = -^--^-l^o. (4) 

n n 

At zero temperature D(T = 0) = Dq has been introduced || to distinguish an ideal conductor with Dq > from 
an insulator with Dq = 0. Our aim here is to analyze the transport behavior at finite temperatures. For orientation, 
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at T > 0, a conductor can develop either to a normal conductor (resistor) with D(T) — but <jq = a(u> — > 0) > 0, or 
remain an ideal conductor characterized by D(T) > 0. An insulator might develop to a normal conductor (conducting 
by thermally activated transport) with D(T) = 0, 0o(T) > 0, remain an ideal insulator with D(T) = 0, oo(T) = or 
even become an ideal conductor with D(T) > 0. 

Below we present results for a(u) for the prototype ID tight-binding model of interacting spinless fermions with 
nearest and next-nearest neighbor interaction. For systems with Hilbert space dimension less than typically 1000 states 
(after implementation of translational symmetry) we calculate o~(u>) directly from Eq.(|^) by finding all eigenstates and 
evaluating current matrix elements; for systems with larger basis dimensions we use a T > Lanczos-based numerical 
technique ||. 

As we are interested mostly in differences in the qualitative behavior of integrable vs. non-integrable systems we can 
restrict our study to high temperatures thus minimizing spurious effects due to the sparse low energy level spectrum 
in finite size systems. It corresponds, in normal conductors, studying systems with mean free paths shorter than the 
lattice size. 

Further, we present the integrated and normalized intensity 

9 T C u or n 

J(w) = D* + \ duj'aiuj')- D* = (5) 

(-T) J 0+ (-T) 

as it contains the relevant information in the conductance and avoids the smoothing procedure of the discrete o~(u>) 
spectra of finite size systems. 

The Hamiltonian we study is given by: 

L L L 

i—1 i—1 i—1 

where Ci{c\) are annihilation (creation) operators of a spinless fermion at site i, rii — c\c\. This Hamiltonian is 
integrable using the Bethe ansatz method for W — [§J and non-integrable for W ^ 0. For W — and V < 2t the 
ground state is metallic, while for V > 2t a charge gap opens and the system is an insulator. 

We study numerically various size systems with periodic boundary conditions and M — L/2 fermions (half- filling). 
The results for L = 8, 12, 16 are obtained by the complete diagonalization of the Hamiltonian, while for L = 20, 24 
the Lanczos method is employed. It should be mentioned that in the latter cases results, e.g. for D* , are subject to 
small statistical error due to finite random sampling Q. 

Metallic state: In Fig. 1 we show the finite temperature conductance for an integrable case. To study the finite size 
dependence of the charge stiffness, we plot in the inset D* as a function of 1/L; the dashed lines indicate a 3 rd order 
polynomial extrapolation based on the L = 8, 12, 16 site systems, suggested by the very good agreement obtained 
with the T = analytical result (square at 1/L — 0) ||. We find that for L — » oo the extrapolated D* ^ 0. At the 
same time cto = o~(lo — > 0) — > as I{uj) seems to approach uj = with zero slope (<t(uj) is the derivative of I(u>))- This 
behavior is reminiscent of a pseudogap. These two results indicate that the integrable system behaves as an ideal 
conductor at T > 0. Moreover we find that the normalized D* approaches a nontrivial finite value in the limit 
T — > oo, depending on V/t and filling, as both D and (— T) are proportional to (3 in this limit. 

In Fig. 2 we show I(oj) and D* for a non-integrable case. Here, as expected for a generic metallic conductor 
(resistor), we find that D* scales to zero, probably exponentially with system size, and oo > as I(u>) approaches 
u> = with a finite slope. These two results imply that the non-integrable system behaves as a normal conductor at 
T>0. 

To further point out the particularity of integrable systems, we investigate the behavior of the conductance on 
approaching the integrable point W = 0. In Fig. 3 we present I{u>) scanning the parameter W . We clearly recognize 
a continuous transfer of the 8— function weight I(u> = 0) = D* at W = to low frequencies, both for W > and 
W < 0. From a calculation of the second frequency moment of the conductance at infinite temperature, we estimate 
the frequency range of a(uj) > proportional to (V — W) 2 + W 2 /2t 2 . Due to remaining finite size effects we are not 
attempting yet to make more quantitative statements about the critical behavior of the low frequency conductance. 

These numerical results on the finite temperature charge stiffness, although not conclusive, strongly suggest a 
qualitative difference between integrable and nonintegrable systems. We can argue about this difference by considering 
the expression for D as a thermal average of the curvatures of levels subject to a fictitious flux. The integrable systems, 
as they are characterized by absence of repulsion between levels, they allow larger fluctuations in the level response 
to a flux and thus plausibly a finite charge stiffness. 

On the other hand in the nonintegrable systems, because of level repulsion, the motion of levels with flux is 
constrained within a characteristic level spacing that is proportional to e~ aL , the density of states, and thus to a 
vanishing charge stiffness with increasing system size. 
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We have also verified that in our integrable system the absence of level repulsion leads to Poisson statistics of the 
level spacings while in the integrable one the level repulsion leads to GOE statistics Jicj . 

As for a o . it is more difficult to ascertain its behavior in the infinite size limit from numerical results in finite 
size systems. For the nonintcgrable systems we find, as expected, a finite value for ao- For the finite size integrable 
systems that we can study, although I(u>) seems to approach u) = with a zero slope, we cannot really exclude a 
finite slope for L — > oo. However, from the physical point of view, even in this case one can expect ideal conductance 
provided the charge stiffness remains finite. It is indicative of a free accelerating system similar in the spirit of a two 
fluid model. 

Insulating state: In Fig. 4 we show I(oj) for the integrable case V = 4t,W = 0. At this value of V > 2t the ground 
state is insulating characterized by D = and a charge gap A ~ t [0. We find that, at finite temperatures, 
D*(T > 0) = I(uj = 0) seems also to decrease exponentially with the system size scaling to zero for L — > oo. This 
precludes a possibility for ideal conductance at T > 0. Furthermore I(u>) seems to approach u) = with zero slope, 
showing a depletion of weight within a low frequency region of order to < t. These are characteristics of an ideal 
insulator, not conducting even at high temperatures T ^> Aq. 

In contrast, as shown in the inset of Fig. 4 (for L = 16), non-integrable systems of roughly the same charge 
gap exhibit a qualitatively different behavior. I(lo) approaches lo — with finite (although small) slope, consistent 
with a small static conductance ao > 0. As expected, conductance here is of a thermally activated character, since 
appreciable ao > appears only at elevated T 3> t. 

To discuss the behavior of the conductance in the insulating regime, it is a good starting point to think about 
the large V/t limit. In this limit the energy spectrum splits in "Hubbard" bands with a fixed number N s of soliton- 
antisoliton (ss) pairs, solitons corresponding to nearest neighbor occupied sites while antisolitons to nearest neighbor 
empty sites. In the V = oo limit solitons and antisolitons are impenetrable. Crossing can only occur by annihilation 
and creation of a pair which corresponds to mixing with other "Hubbard" bands and it has an amplitude of order 
t 2 /V . The states are grouped in characteristic sequences of solitons and antisolitons and the energies of the eigenstates 
are the same as those of a gas of 2A^ S impenetrable particles. 

To analyze the flux <p dependence of the energy we note that the phase associated with the hopping of a soliton in 
a given direction is opposite to the phase picked by an antisoliton; so solitons and antisolitons carry opposite charge. 
It follows that, in this V = oo limit, the flux <fi dependence of the hopping matrix elements can be removed by a 
redefinition of the phase of the states. An nontrivial <fr dependence would have appeared if by successively applying 
the Hamiltonian on a given ss state we could bring it to an equivalent one but with an accumulation of a nonzero 
phase factor. This process, which corresponds to a uniform translation of fermions, is not possible provided we do 
not allow for soliton-antisoliton crossing. 

Therefore, as we have also numerically verified, the width of the "Hubbard" band is of the order of N s t but the 
energy levels are independent of the flux <j> and D is strictly zero. At the same time it is reasonable to argue that 
a static electric field acting on impenetrable particles of opposite charge cannot produce a uniform current and the 
static conductance ao should also be zero. Thus, in this V — oo limit, we find an ideal insulator at any temperature 
independently of the integrability of the system. 

Now, allowing crossing of solitons and antisolitons leads to a D ^ in a finite size system. However our numerical 
results above suggest that D scales again to zero exponentially with the system size. To estimate the dependence 
of the energies in the flux <p due to ss crossing we can argue in terms of ordinary perturbation theory. As above, 
a (j> dependence will appear if by applying the Hamiltonian on a ss state we can bring it to an equivalent one, 
accumulating in the process a nonzero phase factor. This process, that is now possible, necessitates the successive 
creation and annihilation of N s ss pairs and will appear in the ~ N s order in perturbation expansion. As it involves 
N s intermediate states with energy difference V it is an exponentially small contribution. 

Further, taking the point of view that the large V/t limit is the fixed point of the insulating behavior, we can argue 
against the scenario of an insulator at zero temperature developing to an ideal conductor at finite temperatures. 

The next point to discuss is the possibility of an integrable insulating system developing into a normal conductor 
(resistor) at finite temperatures. Unfortunately this point can only be settled after clarifying the exact connection 
between the static conductance, ao, and integrability. However, consistently with our conjecture, we have found 
from independent spectral analysis that whenever the insulating system is integrable, the level spacing distribution is 
Poisson and ao seems to vanish, while when the system is nonintcgrable, the statistics is GOE and ao ^ 0. 

Finally, in order to visualize the above picture of the insulating state, we present in Fig. 5 calculations in the 
large V/t limit. For V = 8t the system is characterized by a much larger charge gap Ao ~ 6t. Indeed we see a 
region of finite conductance in the frequency region < u> < At corresponding to excitations within the first, one 
soliton-antisoliton pair, "Hubbard" band. The weight in this region is increasing exponentially with temperature, a 
sign of activated transport. It is followed by a vanishing conductance up to w ~ 6i, when transitions from the ground 
state to the first Hubbard band start. 

Comments: Our conclusions should hold for other integrable systems as well. Since the anisotropic (and isotropic) 
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spin-1/2 Heisenberg model is equivalent to the integrable (W — 0) model (g), analogous conclusions should apply for 
the spin stiffness and spin diffusion at T > 0. Furthermore, we have numerical evidence (to be presented elsewhere) 
that also the integrable ID Hubbard model, exhibits the same features found for the prototype model @. Namely, 
out of half filling, the system seems to be an ideal conductor, while at half-filling results are consistent with an ideal 
insulator with ao(T > 0) = 0, for any strength of the repulsive interaction. 

We should stress that the above results are only indicative of a relation between integrability and finite temperature 
conductance and our arguments are far from rigorous. Further analytical and numerical studies are necessary to prove 
the validity of our conclusions. However, taking into account present limitations on the size of the systems that can 
be numerically studied and the absence of analytical work, we think that the results presented here are qualitatively 
clear enough to warrant further work on this novel idea. Furthermore, an experimental effort is necessary to observe 
this unusual conductance enhancement. Finally the stability of this effect to deviations from integrability should be 
studied, similarly to the problem of stability of classical soliton systems. 
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FIG. 1. Integrated conductance I(ui) for V = 1.54, W = 0, T = 2t, for L — 12 — 16 (exact diagonalization - full fines) and 
L — 20, 24 (Lanczos method - dotted lines). Inset shows normalized charge stiffness D* vs. 1/L: exact diagonalization (disks), 
T > Lanczos method (circles), analytical result (square) and 3 rd order polynomial extrapolation from L — 8, 12, 16 (dotted 
line) . 

FIG. 2. Integrated conductance for V = 1.5*, W = t,T = 2t, for L = 12 - 20. In the inset InD* vs. L is plotted. 
Notation is as in Fig. 1. 

FIG. 3. I{lo) for L = 16, V = 1.5*, T = 2*, and W/t = -0.5, -0.3, 0., 0.3, 0.5. 
FIG. 4. I(uj) for V = 4*, W = 0, T = At, for L — 12, 16 (full lines with increasing line thickness-exact diagonalization) and 
L — 20 (Lanczos method-dotted line). Inset: for L — 16, T = 4t. 

FIG. 5. for V = 8t, W = 0, T = At, for L = 12, 16 with the notation as in Fig. 4. 
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